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Abstract 



The characteristics of the resonant disturbing function for an asteroid perturbed by a planet 
in circular orbit are discussed. The location of the libration centers and their dep endence with the 
orbital elements of the resonant orbit are analyzed. A proposed numerical method (jGallard o 2006a) 
for computing the strengths of the resonances is revised and applied to the region of the main belt 
of asteroids showing the relevance of several mean motion resonances (MMR) with several planets. 
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Introduction 



There is an interesting diversity of studies on asteroids in mean motion resonances (MMR) 
with Jupiter and on transneptunian objects in MMR with Neptune. However, not all possible 
resonances were analyzed neither all the perturbing planets were considered. 

The Solar System is in fact covered by a innumerable quantity of possible resonances. If we do 
not have a method that adequately weighs the strength of each resonance it is laborious to identify 
which one of the hundreds of possible MMRs that theoretically exist near the semimajor axis of 
the orbit we are studying is the one affecting the body's motion. 

We know by basic celestial mechanics (see also section 2) that the strength of a resonance is 
approximately proportional to the the mass of the planet and to the eccentricity of the resonant orbit 
elevated to the order of the resonance, whenever the eccentricity is not very high. Consequently, 
in general we are not motivated to consider high order resonances neither low mass perturbing 
planets. 

For zero inclination orbits and not very high eccentricities it is possible to compute the widths 
(somehow related t o the strengths) in semima j or axis of the MMRs with the planets as a f unction 
of the eccentricity (jDermott and Murraylll983l : iMorbidelli et q/JI 19951 ; iNesvorny et al.ll2002l ) but no 
simple theoretical method exists to compute the widths in the case of very eccentric and non zero 
inclination orbits. 



Recently iGallardol (|2006al ) presented a numerical method to estimate the strength of an arbi- 
trary mean motion resonant orbit with arbitrary orbital elements assuming a circular orbit for the 
perturbing planet. Based on this principle it is possible to compute the strength of the resonances 
with all the planets form Mercury to Neptune for all ranges of semimajor axis, from the Sun up to 
the limits of the Solar System. That tool allowed the author to identify candidates to be in exotic 
resonances like 6:5 and 1:2 with Venus and 1:2 and 2:5 with Earth. Ve ry recently a n umerous 
population of asteroids in the resonance 1:2 with Mars was also identified (jGallardd 120071 ) . 



We resume here the principles of the method, we analyze some consequences and we apply it 
to the main belt of asteroids putting in evidence the signatures of some MMRs. 



2. The Disturbing Function 

Given a planet of mass rap and radius vector vp in an heliocentric frame and a small body at 
r the equation of motion is given by 

f + fe 2 M -^ = VM (1) 



where R is the disturbing function: 

= k 2 mp[ - 

r p — r | rp 



^pfn-^-T) ( 2 ) 
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In order to construct an analytical theory for the dynamics of a small body we need an 
expression for R. Since Laplace's times classical expressions were constr ucted as series expansion s 
around e = and i = 0° being one of the most recent versi ons the one of Ellis and Murravl (|2000l ). 
These expansions have convergence problems for high e (|Ferraz-Melld I1994 ) and it is necessary 
to take into account several terms for properly account for high inclination o rbits. Planar orbits 
with high eccentricities are best handled by other expansions like Beauge's (jBeaugel 1 19961 ). For 
the general problem a local exp ansion can be constructed in order to study t he motion around a 



small region of the phase space (jFerraz-Mello and Sato 



1989 



Roig et aZ.lll998l ). These asymmetric 



expansions are valid around a small region near the center of the expansion but they can be applied 
to very high eccentricity orbits allowing the calculation of the precise positions of the libration 
centers and the periods of the small amplitude librations. 

Considering a system composed by the Sun, the planet and a small body with orbital elements 
(a, e, i, w, O) the classical expression of the expansion for R is a series of terms of the form: 



R = ^Ccos(</?) (3) 

being C a function of the form 

C = A{a)efe k2 sfs M (4) 

with s = sin(i/2), A(a) being a function of a = (a/ap)^ 1 with positive exponent for resonances 
interior to planet's orbit and negative for exterior ones. The angle ip is defined as 

<P = jiXp + j 2 X + fawp + j4_w + j 5 0p + j 6 fl (5) 

where subscript P denotes the planet. The ji are integers verifying = with j§ + jg being 
always even (D'Alembert rules) and fcj >| |. If (p is a quick varying angle the total effect of the 
corresponding term of R will vanish in the long term evolution. Also if cp is a slow varying angle but 
C is vanishingly small the term again will not have a dynamical effect in the motion of the particle. 
Then if we are interested in a correct description of the long term dynamical evolution of the particle 
we must take into account all the slow varying terms with non negligible coefficients C; these are 
the resonant and secular terms. As ep,e,sp,s are less than 1 we say that the corresponding term 
is of order (ks + + k$ + he). It is possible to simplify the analysis taking a circular orbit for the 
planet with zero inclination, that means ep = sp = 0. With this approximation the number of 
terms involved in ([3]) drop considerably. It is also possible to take into account in R the terms due 
to all planets. 

Under these hypothesis and taking into account D'Alembert rules a g-order resonance \p+q\ : \p\ 
with p and q integers occurs when the general critical angle 

a j = ip + <?)Ap — pX — iq — 2j)w — 2j£2 = a + 2ju (6) 

librates or have a slow time evolution, where 



° = ip + Q)Xp — p\ — qvo 



(7) 
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is the principal critical angle and j is an integer positive or negative. Due to the very slow time 
evolution of the angles (w, O) the librations of Oj occur approximately for 



n p + q 

Tip ~ p 

where the re's are the mean motions. Then, the formula 



(8) 



--(i+-pr i/3 HH 2/3 o) 

ap \p + qJ 

defines the location of the resonances with planet P in semimajor axis. At very low eccentricities 
the time variation of w cannot be ignored and the location of the resonances are shifted (the 
law of structure) respect to Eq. ([9]). The integer p is the degree of the resonance with p < 
for exterior resonances and p > for interior resonances. With this notation the trojans (or co- 
orbitals) correspond to p = —1 and q = 0. The resonant motion is generated when there is a 
strong dependence of M. on a which must be librating or in slow time-evolution. In this case M(cr) 
dominate the time evolution of the small body's orbit. 



3. Properties of the Resonant Disturbing Function 

The limitations imposed by the problems of convergence of the analytical de velopments mot i 



vate the authors to explore the disturbing function numerically. In analogy with ISchubartl ([1968 ), 
in order to explore numerically the function M for a specific resonance defined by a semimajor axis 
given by Eq. ([9]) we can eliminate all short period terms on M computing the mean disturbing 
function 

^ r2iv\p\ 

R {°) = 7m R(Xp,X(Xp,a))dX P (10) 

27r M Jo 

for a given set of fixed values of (e, i, w, O, a) where we have expressed A = A(Ap, a) from Eq. ([7]) 
with a as a fixed parameter and where R(Ap, A) is evaluated numerically from Eq. ([2]) where rp 
and r were expressed as functions of the orbital elements and mean longitudes Ap and A. We repeat 
for a series of values of a between (0°,360°) obtaining a numerical representation of the resonant 
disturbing function R(a). 

The equations of the resonant motion show that the time evolution of the semimajor axis, 
da/dt, is proportional to dR/da then the shape of R(o~) is crucial because it defines the location of 
stable and unstable equilibrium points. For specific values of (e, i, uj) the minima of R(o~) define the 
stable equilibrium points also known as libration centers around which there exist the librations, 
that means oscillations of the critical angle a. The unstable equilibrium points are defined by the 
maxima. 



At low (e, i) the function R(cr) calculated from Eq. (jlOf) is very close to a sinusoid with 
amplitude proportional to e q as one can expect from the classical series expansions. At higher 
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eccentricities the orbit approaches to the planet's orbit, R(o~) start to depart from the sinusoid 
and classical series expansions start to fail with some exceptions like Beauge's expansion. For 
eccentricities greater than the collision eccentricity e c : 



I _ , p + q\ 2 / 3 
P 



(11) 



the orbit can intersect the planet's orbit and for low inclination orbits two peaks start to appear 
around the point where R(o~) has its maximum for e < e c . If the two peaks can be distinguished 
then a stable equilibrium point appears between them. For high inclination orbits the intersection 
between orbits is less probable and soft maxima can appear instead of the peaks. For small inclina- 
tion orbits the libration centers are almost independent of and can be classified in only three 
different classes as showed at Table [TJ 

Resonances of the type l:n including trojans (that means 1:1) exhibit a similar general behavior. 
For e < e a , where e a is certain value verifying e a < e c , th ere is a librati on center at a = 180° and 



for e > e a there appear the asymmetric libration centers (lBeaugelll994l ) with locations depending 
not only on e but also on (i,to). For trojans we have e a = e c = so the low eccentricity librations 
around a = 180° do not exist. For e > e a horseshoe (HS) trajectories wrapping the asymmetric 
librations can exist. These HS trajectories are of the same nature of the horseshoe trajectories 
in the case of trojans and are only possible for this kind of resonances. In HS trajectories a is 
oscillating with high amplitude around 180°. For e > e c at low inclinations they appear two peaks 
(unstable equilibrium points) and a stable libration center at a = 0°. T his last equilibrium point 



is associated with the known quasi- satellites (QS) of the 1:1 resonances (jWiegert et al\ 12000 ) 



For low inclination orbits all odd order interior resonances show librations around a = 0°, and 
for e > e c it appears another libration point at a = 180°. Conversely, all interior resonances of order 
even and all exterior resonances except resonances of type l:n show librations around a = 180°, 
and for e > e c it appears another libration point at a = 0°. 

For high inclination orbits the geometry of the encounters is strongly modified affecting the 
peaks of R(cr) and its shape becomes completely different to the low inclination case becoming the 
libration centers strongly dependent on oj. Figures HJ El and [3] illustrate the dependence of R(o~) 
and its libration centers with (e,i,u) for the case of the exterior resonance 1:4 with Neptune. In 
all figures we have taken such units that k 2 mj up = 1. 



4. The Strength of the Resonances 

For a given resonant orbi t defined by para meters (e, i, w, O) with a given planet the disturbing 
function R(o~) is determined. iGallardol (|2006al ) defined the strength function SR as: 



SR(e,i,uj) =< R > -R Tl 



(12) 
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being < R > the mean value of R(cr) with respect to a and R m i n the minimum value of R(a). 
This definition is in agreement with the coefficients of the resonant terms of the expansion of the 
disturbing function for low (e, i) orbits because for this case R{a) is a sinusoid with an amplitude 
given by < R > —R m in- Then, for low eccentricity and low inclination orbits SR should follow 
the function e q . For high (e, i) resonant orbits the strength cannot be calculated by analytical 
developments and departures form the low eccentricity regime is the rule. Figures [H [5] and [6] 
illustrate the dependence of SR with (e, i, lo) for the case of the exterior resonance 1:4 with Neptune. 

When SR ~ we have dR/da ~ for all values of a and then da/dt will not be dominated 
by resonant terms but by other terms that will generate some time evolution of the semimajor 
axis and consequently the resonance will be broken, then the resonance will not be dynamically 
significant or strong. On the contrary, a high value of SR implies a strong dependence of R on 
a and the resonant disturbing function R(cr) will dominate the motion forcing the semimajor axis 
to evolve oscillating around the stable equilibrium points or to evolve escaping from the unstable 
equilibrium points. 



Gallardd (|2006al ) analyzed the shape of SR(e,i,uj) for several resonances and found that all 



them can be roughly classified in two classes (q < 1 and q > 2) according to the response of SR 
to the variation of the inclination which is more important for resonances of order 2 or greater. 
It is possible to understand why the inclination is an important factor for resonances of order 2 
or greater. Analytical developments of R(cr) in powers of (e,i) show that for a q-o rder resonance 
the lowest o rder resonant terms are of order q in (e,i) ([Murray and Dermottlll999l ). In particular 
for trojans (jMoraislll999l ) and first order resonances the lowest order terms are independent of i. 
But, for resonances of order q > 2 the lowest order resonant terms have a dep endence with i that 

In consequence is 



make some contribution to R(<r) for high inclination orbits (jGallardo 



2006 



epe 

3). 



natural that for resonances of order 2 or greater the resonances are stronger for high inclination 
orbits because the resonant terms depending on i will show up. On the contrary we cannot expect 
such behavior for resonances of order 1 or because the resonant terms depending on i have lower 
relevance. 



5. Applications: Identification of the Resonant Signatures 



Following the numerical procedure we have described, iGallardd (|2006al ) presented an atlas of 
MMRs for the Solar System and identified several objects in unusual resonances with the terrestrial 
and with the jovian planets. In particular it was found there exist asteroids evolving in the exterior 
resonances 1:2 and 2:5 with the Earth. 



We present here an atlas of resonances for the main belt of asteroids reworked from iGallardc 



(|2007l ) and we compare it with the distribution of known asteroids (Fig. [7]). The atlas was cal- 
culated from Eq. (|12p considering e = 0.3, i = 10° and to = 60°. The histogram of asteroids 
was elaborated using the osculating orbital elements of ~ 370000 asteroids taken from ASTORB 
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database (ftp://ftp.lowell.edu/pub/elgb/astorb.html) and grouped in bins of 0.001 AU. It is pos- 



sible to identify the signatures of some well known resonances with Jupiter. In general, a secular 
evolution inside these resonances drives the eccentricity to values that a collision with the Sun or a 
close encounter with Mars or Earth remove the asteroid from the resonance so a gap is generated. 

The main belt of asteroids is limited at its extremes by the resonances 4:1 (more precisely a 
secular resonance is responsible for this border) and 2:1 with Jupiter. The depletion effects due to 
resonances 3:1, 8:3, 5:2, 7:3, 9:4 and also 11:5 with Jupiter are evident. But we can also focus in 
some signatures not so evident. It is possible to identify some excess of asteroids in the location of 
certain resonances. In particular, in the histogram there is an excess of around 30% of asteroids at 
a ~ 2.419 AU, exactly where the resonance 1:2 with Mars is located. Note also that the resonance is 
isolated and consequently not perturbed by others, that means, it should dominat e in that r egion . 



That population inside the resonance was confirmed via numerical integrations by lGallardd ((2007) 
constituting around a thousand of known asteroids evolving in the resonance. This is the first 
numerous population that we have knowledge captured in a MMR with a terrestrial planet. 



Figure [7] is also showing that resonances 2:5E already studied by iGallardol (|2006al ) and 3:4M, 
3:8E, 2:3M and 1:3E are strong enough to be considered possible reservoirs of asteroids, at least 
temporarily. This last one probably is strongly perturbed by 4:1 J but 2:3M is relatively isolated 
and 3:8E is at the middle of a considerable population of asteroids, in consequence they should be 
populated. 

As a closing comment, it is evident that analytical theories plus numerical procedures give us 
at present a quite complete understanding of the MMRs and several features of the distribution of 
asteroids' populations in the Solar System can be understood in this context. 
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Resonance Type 




new do for e > e c 


exterior l:n 


180° or asymmetric 


0° 


exterior others 

and 
interior q even 


180° 


0° 


interior q odd 


0° 


180° 



Table 1: Stable equilibrium points for low inclination resonant orbits. Resonances of type l:n except 
trojans have an equilibrium point at a = 180° for low eccentricity orbits. For higher eccentricities 
this point bifurcates in the two asymmetric points. For eccentricities greater than the collision 
eccentricity, e c , it appear another new equilibrium point in all resonances. For high inclination 
orbits this scheme is strongly modified and the argument of the perihelion becomes relevant for the 
location of the equilibrium points (see Fig. [3]). 




Fig. 1. — Resonant disturbing function for the resonance 1:4 with Neptune calculated from Eq. 
(|10p for three different values of the eccentricity. In all cases i = 30° and uj = 60°. The minima 
define the location of the stable equilibrium points. Shallow minima at e = 0.1 implies low stability. 
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Fig. 2. — Resonant disturbing function for the resonance 1:4 with Neptune calculated from Eq. ()10p 
for three different values of the inclination. In all cases e = 0.3 and u = 60°. The minima define 
the location of the stable equilibrium points. Shallow minima at i = 178° implies low stability. 
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a 

Fig. 3. — Resonant disturbing function for the resonance 1:4 with Neptune calculated from Eq. 
()10p for two different values of the argument of the perihelion. In all cases e = 0.3 and i = 30°. 
The minima define the location of the stable equilibrium points. 
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Fig. 4. — Dependence of the strength SR(e) with the eccentricity for the resonance 1:4 with Neptune 
calculated from Eq. (jl2h . At low eccentricity regime SR oc e q with q = 3 in this resonance. 
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Fig. 5. — Dependence of the strength SR(i) with the inclination for the resonance 1:4 with Neptune 
calculated from Eq. fll2|). The strength is in general greater for high inclination (but direct) orbits. 
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Fig. 6. — Dependence of the strength SR(lo) with the argument of the perihelion for the resonance 
1:4 with Neptune calculated from Eq. (|12p . The argument of the perihelion affects the location of 
the equilibrium points but its influence in the strength is the less important. 
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Fig. 7. — Atlas of the strongest MMRs with all the planets in the region of the main belt of 
asteroids where strengths were calculated from Eq. (|12j) assuming e = 0.3, i = 10° and lo = 60°. 
Superimposed is showed an histogram of semimajor axes constructed with bins of 0.001 AU. The 
peak a t resonance 1:2 with Mars is clearly distinguished. This figure was recomposed from lGallardo 



